Nonequilibrium transitions in complex networks: a model of social interaction 
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We analyze the non-equilibrium order-disorder transition of Axelrod's model of social interaction 
in several complex networks. In a small world network, we find a transition between an ordered 
homogeneous state and a disordered state. The transition point is shifted by the degree of spatial 
disorder of the underlying network, the network disorder favoring ordered configurations. In random 
scale-free networks the transition is only observed for finite size systems, showing system size scaling, 
while in the thermodynamic limit only ordered configurations are always obtained. Thus in the 
thermodynamic limit the transition disappears. However, in structured scale-free networks, the 
phase transition between an ordered and a disordered phase is restored. 

PACS numbers: 87.23.Ge, 05.50.+q 



I. INTRODUCTION 



Lattice models are a powerful basic instrument in the 
study of phase transitions in equilibrium Statistical Me- 
chanics, as well as in non-equilibrium systems Q]. Tradi- 
tionally, equilibrium phase transitions have been studied 
in regular lattices, with the critical temperature being 
a non-universal quantity that depends on the particular 
lattice under consideration, while critical exponents and 
some amplitude ratios are universal quantities depend- 
ing only on spatial dimension and some symmetries of 
the order parameter. The detailed structure of the reg- 
ular network connections is, in most cases, irrelevant in 
the sense of the renormalization group. However, recent 
research in the structure and topology of complex net- 
works H has shown that social interactions and, more 
generally, biological and technological networks, are far 
from being regular as well as being also far from a random 
network or from a mean-field network linking all-to-all. 
This has triggered the study of standard models of Statis- 
tical Mechanics in these complex networks. In particular, 
recent studies of the Ising model in the so-called small- 
world IJ] and the scale-free networks have shown that 
the behavior of the model differs from that observed in a 
regular network. 

A small world network Q is generated by rewiring 
with a probability p the links of a regular lattice by long- 
distance random links. The presence of a small fraction 
of "short cuts" connecting otherwise distant points, dras- 
tically reduces the average shortest distance between any 
pair of nodes in network, keeping the clustering high. The 
small-world networks generated by rewiring links have 
degree distributions with exponential tails. In contrast, 
scale-free networks are characterized by a fat-tailed 



'Electronic address: 


klemm@nbi.dk 




tTIRT,: |ht.tp://HW. 


Tnedea.uib.es/ 


'"victor 


victor@imcdoa.uib.cE 






-f Electronic address: 
§ Electronic address: 


•aul@imedea.uib.eE 




TLaxi@imedea.uib.e 





Electronic address: 



(power law) degree distribution. They are usually mod- 
elled by growing networks and preferential attachment of 
links. 

The Ising model in small-world topologies shows a 
change of behavior from the regular case to the mean 
field characteristics. In Ref. it is shown analytically 
that for a small world lattice, obtained from rewiring 
with probability p the links of a 1-d ring lattice with 2k 
nearest neighbors interactions, the crossover temperature 
to the mean field critical behavior varies for p << 1 as 
Tco{p) oc —k{k -\- l)/ln(p) whereas the critical temper- 
ature scales as Tc{p) oc —2k/\n{p), so that a ferromag- 
netic ordered phase exists for any finite value of p. The 
crossover to mean field behavior in small world ring lat- 
tices has been further discussed in Refs. 0, ^, whereas 
numerical results in 2-d and 3-d lattices have been re- 
ported in j^. Interestingly, if directed links are consid- 
ered, not only the critical temperature changes but the 
nature of the transition also switches from second to first 
order Q. 

A much different behavior is observed in scale-free net- 
works. This can be related to the influence of the pres- 
ence of so-called hubs, i.e. units whose degree is much 
larger than average. This is well illustrated by the behav- 
ior of the Ising model in scale-free networks with degree 
distribution P(fc) oc fc-^, 7 > 1. The results of Refs. 
and [p^ show that equilibrium systems exhibit a phase 
diagram that is qualitatively different from the mean-field 
case. In particular, the Ising model in a random scale- 
free networks shows an infinite critical temperature in the 
thermodynamic limit of an infinite number of nodes. In 
fact it has been developed analytically a mean field the- 
ory connecting the exponent of the degree distribution 
and the critical behavior of the Ising model ^ ^ . 

In this paper we address the question of the role played 
by the topology of complex networks in non-equilibrium 
transitions of models in which there is interaction be- 
tween the variables associated with the nodes connected 
by links in the network. This is a natural next step be- 
yond the analysis of equilibrium, Ising-type models in 
these complex networks. Simple non-equilibrium models 
closely related to percolation have been already consid- 



ered |T^, |l^. Here, and given the social motiva- 
tion and relevance of these complex networks, we have 
chosen to analyze the modelproposed by Axelrod for 
the dissemination of culture |19|] . The spreading process 
in this model cannot be reduced to a percolation pro- 
cess. The model rather describes a competition between 
dominance and spatial coexistence of different states in 
a non-equilibrium dynamics of coupled Potts-type mod- 
els. The model was originally considered by Axelrod in 
a square lattice. The Statistical Mechanics analysis of 
the model in this regular two-dimensional network identi- 
fies a nonequilibrium order-disorder phase transition [ pO[ . 
However, it is interesting to notice that, in his original pa- 
per, Axelrod already discussed the relevance of the topol- 
ogy, speculating that "With random long-distance inter- 
actions, the heterogeneity sustained by local interaction 
cannot be sustained." In particular we consider here this 
question. 

In the next Section we introduce the original model 
proposed in Ref. and summarize briefly the main re- 
sults in regular 2-d networks. The model in small-world 
and scale- free networks is presented in Section |l| and IV 
respectively. The non-equilibrium transition is shown to 
disappear in the thermodynamic limit of a scale free net- 
work. We then consider in Section ^ a structured scale 
free network ||2l[ which shares characteristics of small- 
world and scale-free networks. A non-equilibrium transi- 
tion is shown to persist for large systems in this network. 



Our conclusions are summarized in Section Vl 



II. THE MODEL 

The model we study is defined by considering 
N agents as the sites of a network. The state of 
agent i is a vector of F components (cultural features) 
((Til, (Ti2, ■ ■ ■ , Cif)- Each aif can take any of the q inte- 
ger values (cultural traits) 1 , . . . , g, initially assigned in- 
dependently and with equal probability 1/q. The time- 
discrete dynamics is defined as iterating the following 
steps: 

1. Select at random a pair of sites of the network con- 
nected by a bond (i, j). 

2. Calculate the overlap (number of shared features) 

3. If < ^(i, j) < F, the bond is said to be active and 
sites i and j interact with probability l{i,j)/F. In 
case of interaction, choose g randomly such that 
'^ig 7^ '^■jg ^^^d set = cr.jg. 

In any finite network the dynamics settles into an 
absorbing state, characterized by the absence of active 
bonds. Obviously all the completely homogeneous 
configurations are absorbing. Homogeneous means here 
that all the sites have the same value of the cultural trait 
for each cultural feature. Inhomogeneous states consist- 
ing of two or more homogeneous domains interconnected 
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FIG. 1: The order parameter Smax/N as a function of q in 
regular lattices with A'' = lOO'^ nodes for F = 10. For each 
value of q the outcome of 32 independent realizations is shown. 
The transition occurs for g ~ 55 (see Fig.B). 



by bonds with zero overlap are absorbing as well. A 
domain is a set of contiguous sites with identical cul- 
tural traits. We define an order parameter in this system 
po| , ^ as the relative size of the largest homogeneous 
domain Smax/N, being N the number of sites in the net- 
work. 

Previous results have been obtained in square lattices 
with nearest neighbor interaction. A variation of the 
model with initial distribution of traits according to a 
Poisson rather than a uniform distribution shows a non- 
equilibrium order-disorder phase transition where the 
number of traits q plays the role of a control parame- 
ter . The system reaches ordered absorbing states for 
q < qc {Smax = 0{N)) and disordered states for q > qc 
{Smax ^ N). The same type of phase transition occurs 
in the original model with a uniform initial distribution 
of traits p3[ . 

When comparing the effect of complex networks in this 
phase transition with the equilibrium Ising transition one 
should notice several conceptual differences. First, this is 
a sort of zero-temperature transition in which ordered or 
disordered states exists with no reference to thermal fluc- 
tuations. In fact, the effect of small noise in this system 
is essential, revealing the presence of metastable states 
and changing the phase diagram in a nontrivial way [ p3| . 
A second related point is that the control parameter of 
the transition q, is here not a collective property of the 
system as temperature, but rather an ingredient of the 
definition of the system itself. In a way, the transition 
occurs going from one system to another as q is changed. 
Finally, and in reference to critical properties and expo- 
nents, we note that the transition (except for F=2) is of 
first-order type: In figure |l| we plot the final values for the 
order parameter, obtained for 32 different realizations of 
the dynamics. Notice that for q < 50 all the systems 
end up in a homogenous state that basically fills up the 
entire lattice {Smax/N « 1) whereas for g > 60 the max- 
imum homogenous regions obtained are very small. This 
is the order-disorder phase transition discussed before. 
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FIG. 2: Distribution of the order parameter at g = 55, F = 10 
in a square lattice of size A'' = 100^ . The distribution is based 
on 100 independent realizations. 

For 50 < g < 60 we observe bistability in the sense that 
the system settles around any of two mean values for the 
order parameter. This bistability, which is usually asso- 
ciated with first order phase transitions, is clearly made 
explicit in the corresponding histogram shown in figure |^ 
where the two preferred values appear as maxima of the 
histogram. The transition point corresponds to q = qc 
for which these two values are equally probable. 



III. SMALL WORLD NETWORK 

Social networks are far from being regular or com- 
pletely random. However they also share some features 
with them. On the one hand, social networks are known 
to be small ]25| , i.e. any pair of nodes in the network can 
be connected following a number of links much smaller 
than the size of the network. This is also observed in ran- 
dom networks, where the average shortest distance be- 
tween pair of nodes (the so called path length i) increases 
logarithmically with the size of the network £ ^ \nN, 
while in regular lattice in d-dimensions £ '--^ iV^/'^. On 
the other hand, social networks are also known to form 
cliques i.e. groups of nodes highly connected be- 
tween them. "Cliquishness" can be characterized by the 
so-called clustering coefficient C, which is defined as the 
relative number of closed triangles in the network. Regu- 
lar lattices can show large clustering while in random net- 
works C ^ N^^. High clustering and short path length 
define a small world network. 

The first model encornpassing the small world effect 
was introduced in Ref. proposing an algorithm that 
interpolates between a random and a regular lattice. 
First one generates a two-dimensional regular lattice with 
bonds between nearest neighbors and open boundary 
conditions. Then for each bond (ij), with probability 
p detach the bond from node j and attach it to a node I 
instead. Node / is chosen at random with the restriction 
that duplicate and self-connections are excluded. The pa- 
rameter p interpolates between the original regular lattice 
(p = 0, no rewiring) and a network very similar to a ran- 




FIG. 3: The average order parameter {Smax}/N as a function 
of q for three different values of the small world parameter 
p. System sizes are N = 500^ (squares) and N = 1000^ 
(diamonds), number of features F — 10. Each plotted value 
is an average over 100 runs with independent rewiring {p > 0) 
and independent initial conditions. 

dom graph {p = 1). Thus in the limiting case p = we 
have a network with high clustering but also large path 
length; in the limit p = 1 we have networks with small 
path length but also small clustering. For intermediate 
values of p the algorithm generates networks with high 
clustering and small path length. 

We now study the behavior of Axelrod's model in de- 
pendence of p. A small world network is used from the 
beginning of each simulation run. Figure ^ shows the 
dependence of the order parameter on q, for three dif- 
ferent values of p. For any fixed value p > we find 
a nonequilibrium phase transition which becomes sharp 
and well defined as the system size increases. There is a 
critical value qc of the control parameter which separates 
the ordered and the disordered state, just as in regular 
lattices. However, qc increases with the amount of spa- 
tial disorder. This is clearly shown in the (p, 5)-phase 
diagram. Fig. ^. The filled area above the {p,qc{p))- 
curve represents the disordered states, the area below the 
curve represents the ordered states. Consequently, for 
values q < qdp) the outcome of the dynamics is always 
complete order, whereas for q > qdp) only disordered 
frozen states are encountered. The density p of rewired 
bonds determines the nature of these frozen states, but 
for q < qdp = 1) the system orders by increasing p, that 
is, the number of long distance links. We find a depen- 
dence qdp) — qdp = 0) cx p" with a = 0.4 obtained from 
a best fit. This result is displayed in the inset of Fig. 
^. Therefore, we find the same qualitative result that for 
the equilibrium Ising model, in the sense that the small 
world connectivity favors ordered states. 

The robustness of the phase diagram is shown by per- 
forming a different dynamical scenario. First, a run of the 
dynamics in a regular lattice is performed. Only after an 
absorbing configuration has been reached the lattice is 
rewired according to the above rewiring procedure with 




FIG. 5: The average order parameter {Smax)/N in random 
scale-free networks for F = 10. Averages are taken over 1000 
independent realizations. Different curves are for different 
system sizes: 1000 (circles), 2000 (squares), 5000 (diamonds), 
and 10000 (triangles). 



FIG. 4: Phase diagram for the Axelrod model in a small world 
network. The curve separates parameter values (p,q) which 
produce a disordered state (shaded area) from those with or- 
dered outcome (white area). For a given p the plotted value 
qc is the one for which the value of the order parameter is 
closest to the, somewhat arbitrary but small, value 0.1 for 
system size A'' = 500'^ and _F = 10. Inset: After subtraction 
of a bias qdp = 0) = 57, qdp) follows a power law oc p" "^^ 
(dashed line). 



the parameter p. After the rewiring, the configuration 
is not necessarily frozen because the rewiring can intro- 
duce active bonds connecting compatible cultures that 
have been disconnected before. Starting the dynamics 
again, the system may relax to a different absorbing con- 
figuration, which in general is more ordered than the con- 
figuration reached before the rewiring. After this second 
phase of relaxation, the order parameter is measured in 
the absorbing state. We find that the results of this alter- 
native scenario (see Fig. ^ are in good agreement with 
the ones of the above original scenario, starting with a 
small world network in the initial condition. 



IV. SCALE-FREE NETWORKS 

One important ingredient missing in the small world 
networks considered so far is that the degree distribution 
does not show a fat tail. Although it is not clear whether 
social networks present a power law distribution of de- 
gree, the evidence indicates that they_are ubiquitous in 
biological and artificial networks |27, 



28 . Scale-free net- 



works are characterized by a power law tail in the degree 
distribution of the form P{k) oc fc~^ where the exponent 
7 lies in the range between 2 and 3. Two ingredients 



have been shown to be sufficient to generate such feature: 
growing number of nodes and preferential attachment of 
links. 

The well established Barabasi-Albert model is based 
in these two mechanisms At each time step a new 
node is added to the network and attaches m links to an 
existing node with degree k with probability n(fc) cx k. 
This algorithm generates networks whose degree distribu- 
tion follows P{k) = 2m^k~^, the path length increases as 
£ ^ \nN, and the clustering decreases as C ~ (In iV)^/iV. 
We have studied the dynamics of Axelrod's model for the 
diffusion of culture in scale-free networks following this 
algorithm. In Fig. ^ we show our results for the order 
parameter for different system sizes. For a given size N 
we find a transition at qc{N). We can define the critical 
value qdN) as the value where the standard deviation of 
the distribution of Sma.x/N reaches the maximum value. 
We find that qdN) - ArO-39. Using this result we observe 
data collapse with a rescaling qN~^ , see Fig. ^. The best 
result is obtained for (3 = 0.39 in excellent agreement 
with the scaling obtained previously. This indicates that 
in the thermodynamic limit the transition disappears and 
the ordered monocultural state establishes in the system. 
This behavior is similar to the Ising model in regular and 
scale-free networks: While in a two-dimensional lattice 
the Ising model displays a phase transition at a finite 
critical temperature, in random scale-free networks an 
effective transition is observed for finite systems where 
the effective critical temperature diverges logarithmically 
with system size. This can be explained by the role of 
the hubs (nodes with a large number of links) in these 
networks. They are highly instrumental in establishing 
ferromagnetic order in the system. The same prominent 
role is played by the hubs in the case of the dissemina- 
tion of culture. The hubs help the spreading of cultural 
traits as can be inferred from the observed dependence 



0.8 
§ 0.6 - 
0.4 
0.2 







o o I 



o 



10 



qN 



20 

0.39 



30 



40 



FIG. 6: Rescaled plot of the data shown in Fig. ^for different 
system sizes. 
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FIG. 8: Distribution of the order parameter in scale-free net- 
works of size A'^ — 5000 for q — 360, and 1000 realizations. 




FIG. 7; Order parameter in random scale-free networks of 
size N — 5000. For each value of g (i^ = 10) the outcomes of 
100 independent realizations are shown. 



with system size. Note, however that the effective tran- 
sition of Axelrods model for a finite system in a scale 
free network displays the characteristics of a first order 
transition: We show in Figs. ^ and ^ the same type of 
behavior observed in Figs. |^ and |^ for the regular net- 
work. For a range of values of q around qc a realization 
ends either in an ordered monocultural state or in a dis- 
ordered frozen configuration, the two preferred values of 
the order parameter clearly seen in the histogram of Fig. 



V. STRUCTURED SCALE-FREE NETWORKS 

The scale- free networks considered in the previous sec- 
tion underestimate the clustering observed in real net- 
works |2^. A question that merits to be addressed is 
if this low clustering coefficient is responsible for the ab- 
sence of the phase transition in the thermodynamic limit. 




FIG. 9: The average order parameter {Smax)/N as a func- 
tion of q for F = 10 in structured scale-free networks. The 
networks contained A'^ — 1000 (circles), A'^ = 2000 (squares), 
N = 5000 (diamonds) , and TV = 10000 (triangles) nodes with 
F = 10 features. Each data point is an average over 32 inde- 
pendent realizations. 



In order to reproduce a high clustering along with a scale- 
free distribution of the degree, we employ the networks 

Again 



generated by the algorithm proposed in Ref. |21 



at each time step we add a new node to the network. The 
node is attached to the m active nodes in the network. 
The new node becomes active and one of the m + 1 ac- 
tive nodes is deactivated with probability n(fc) oc k~^. 
Starting from m fully inter-connected active nodes, this 
algorithm generates scale-free networks with a cluster- 
ing coefficient C w 5/6 independent of system size. In 
the following we call these networks structured scale-free 
networks because of the large clustering coefficient, the 
strong neg ative correlation between degrees of connected 
nodes |p^ , and the modular structure |2|]. These prop- 
erties are not found in the random scale-free networks of 
the previous section. 

As shown in Fig. p|, in the structured scale-free net- 
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FIG. 10: Relation between the size ol the largest cultural 
cluster and the largest degree in the network for g = 20 (tri- 
angles), q = 100 (squares) and q = 500 (circles). Each data 
point is the outcome of one realization run in a structured 
scale-free network of size A'^ = 1000. Solid lines are running 
overages over 10 adjacent data points for q = 20, 100, 500 (top 
to bottom). For each value of q, 100 independent networks 
and initial conditions were generated. 

works the model displays a behavior different to what we 
observed in random scale-free networks in the previous 
section. For g ^ 10 the system settles into an ordered 
state. For increasing values of q, the order parameter 
undergoes a decay whose slope grows with system size. 
This suggests a phase transition at qc ~ 10, in contrast to 
the absence of a transition point found for the randomly 
wired scale- free networks in the thermodynamic limit. As 
on large scales the structured scale-free networks have 
one-dimensional topology it seems natural that this 
transition at (?c ~ ^ coincides with the behavior of the 
model found in one-dimensional regular lattices p^ . 

At difference with the regular lattices, in the struc- 
tured scale-free networks for q > qc the order parameter 
does not tend to zero. It reaches a finite plateau value, 
indicating partial ordering of the system. Only for val- 
ues q^ qc the order parameter drops below the plateau 
value and tends to zero. This behavior may be under- 
stood by relating the size Smax of the largest cultural 
cluster with the largest degree kmax present in the net- 
work, as shown in Fig. nOl In the intermediate range 
050 < q < 200, where the plateau of the order parameter 
is observed, we find Smax ~ kmax for almost all realiza- 
tions. This suggests that the largest hub, the node with 



the largest degree, and its neighbors order such that they 
form the largest cluster in the absorbing state. As q is 
reduced and its value approaches qc from above, the or- 
dering goes beyond the largest hub and an increasing 
part of the network forms the largest cluster. On the 
other hand, for large values q > 200, the neighborhood 
of the largest hub does no longer reach complete order- 
ing and Srnax < kmax- The q value for the onset of the 
decay of Smax below kmax is expected to be dependent 
on system size: as Fig. ^ shows, for increasing system 
size, the plateau in the order parameter extends to larger 
values of q. With inrceasing system size, the value of 
the plateau is expected to decrease as kmax/N — N^~^ 
with (3 — — 1)^^ , where 7 is the exponent of the de- 
gree distribution. These results suggest that in the limit 
— > 00, the dynamics of the social interaction model in 
structured networks experiences a transition similar the 
one observed in a one-dimensional lattice. 



VI. CONCLUSIONS 

We have found that the nonequilibrium transition be- 
tween order and disorder that exists in a regular d=:2 net- 
work for Axelrods model of cultural influence |l^ is mod- 
ified by underlying complex networks with similar qual- 
itative features that an equilibrium thermal Ising-type 
transition. We have shown that the transition pertains 
also in the presence of random long-distance connections: 
with increasing density of long-distance connections in a 
small world network, the critical point qdp) increases. 
Therefore the small world connectivity favors cultural 
globalization as described by the ordered state. The value 
of qc reaches a maximum for the random network ob- 
tained with a p ~ 1 probability of rewiring in the small 
word network construction. A transition from disorder to 
order is obtained increasing p for a fixed value of the con- 
trol parameter q < qc{p = 1). We have also found that, 
for a fixed finite system size, the scale free connectivity 
is more efficient than the limiting random connectivity of 
the small world network in promoting the ordered state 
of cultural globalization. However, there is a system size 
scaling in the transition observed for a free scale network, 
so that the transition disappears in the thermodynamic 
limit: In the presence of scale-free interactions the order 
state prevails due to the presence of hubs. The consider- 
ation of structured scale free-networks restores the order- 
disorder transitions in spite of the hubs, but the value of 
the order parameter for the disordered state reveals the 
existence of ordered clusters. 
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